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ABSTRACT 


Radar  clutter  spectra  having  sufficient  dynamic  range  and 
resolution  to  permit  shape  estimation  are  reviewed.  All  are  found  to 
be  distinctly  non--Gaussina.  A  clutter  scatterer  motion  model  that  in¬ 
cludes  correlation  of  velocities  is  defined  and  supported  by  empirical 
data  for  the  case  of  wind  fields.  The  Doppler  spectra  of  scatterers 
with  correlated  motions  are  computed.  Improvement  factors  for  binominal- 
weighted  moving-target  indicators  (MTI)  are  computed  for  several  analytic 
spectrum  functions  that  bound  various  observed  features  of  empirical 
spectra.  , 
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1.0 


INTRODUCTION  AND  SUMMARY 


1.1  INTRODUCTION 


In  Barlow's  1949  tutorial  paper  on  Doppler  radar  [l]  he  introduced  the 
notion  of  clutter  spectrum.  He  stated  that  the  power  spectra  of  land,  sea, 
rain,  and  chaff  were  all  approximately  Gaussian, 

S(f )  =  exp  {-  a  f 2 / f Q  }  ,  (1) 

and  he  gave  typical  values  of  the  eoustant  a  for  each.  Since  that  time  , 
almost  all  analysis  of  radars  and  most  interpretation  of  experimental 
observations  of  clutter  proceed  under  the  assumption  of  a  Gaussian  shaped 
spectrum.  The  fluctuation  spectra  produced  in  radar  signals  by  antenna 
scanning  or  by  linear  shear  are  reasonably  expected  to  be  approximately 
Gaussian.  But  we  are  talking  here  about  the  intrinsic  motion  of  the 
scatterers  which  is  the  limiting  spectral  spreading  mechanism  for  noiseless 
non-scanning  radar . 

Since  1950  almost  all  papers  that  have  'Doppler  Spectrum'  in  their 
titles  report  only  the  standard  deviation  of  the  spectrum  width,  and  occasion¬ 
ally  a  few  curves  spanning  10  dB  or  so  of  dynamic  range.  With  so  shallow 
a  range,  no  shape  determination  can  be  made.  But  whenever  spectra  were 
analyzed  to  reveal  spectral  fall  off  to  20  dB  or  more,  the  spectral  shapes 
have  been  distinctly  non-Gaussian.  Rather,  a  shape  more  simply  exponential 
is  suggested: 

(f )  =  exp  {-b 1 f | }  (2) 

This  shape  cannot  be  explained  by  physical  models  of  scattever  motion  in¬ 
volving  either  linear  shear  or  uniform  and  independent  Gaussian  random 
motion  of  scatterers.  Some  other  mechanism  of  correlated  motion  must 
dominate.  Note  that  we  expect  that  motion  must  dominate  in  some  cases, 
rather  than  amplitude  modulation  of  the  scatterers'  cross  sections,  so  a 
motion  based  model  is  first  most  attractive  for  attention. 
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It  is  t'ue  purpose  here  Co  explore  a  scatterer  motion  model  that 
predicts  a  spectrum  more  like  observed  clutter  spectra  than  the  traditional 
Gaussian  assumption. 


The  motion  model  chosen  for  entry  is  a  Wiener  process,  with  the 
radial  velocity  as  the  process  variable  and  the  transverse  position  in  the 
radar  beams  as  the  parameter. 

L. 2  Summary 

Reported  empirical  radar  clutter  spectra  having  dynamic  range  of 
20  dB  or  greater  are  reviewed.  It  is  found  that  none  of  these  spectra  have 
shapes  that  are  reasonably  interpreted  as  Gaussian.  Rather,  they  all  have 
two  features  not  characteristic  of  the  classical  Gaussian  shape.  Their 
shape  near  zero  modulation  frequency  shows  a  peak,  rather  than  a  rounded 
nose,  and  the  falloff  of  spectral  density  at  higher  modulation  frequencies 
is  much  slower  than  Gaussian.  For  the  first  20  uB  or  so  of  density,  the 
shapes  arc  approximately  exponential  in  modulation  frequency. 

A  motion  model  is  postulated  to  account  for  the  observed  spectra, 
Clutter  scatterers  are  assumed  moving  such  that  their  radial  velocities  are 
correlated.  The  correlation  is  introduced  by  assuming  that  the  differences 
in  radial  velocities  of  scatterers  separated  by  a  given  distance  is  a  Weiner 
process  in  the  separation  variable.  Spectra  of  this  model  are  computed  and 
compared  with  the  data  base  for  both  coherent  and  noncoherent  cases.  The 
spectra  computed  from  the  model  have  the  same  sharp  peak  as  the  data,  but  the 
match  to  the  shape  in  the  high-frequency  region  is  indeterminate. 

Several  analytic  functions  having  features  bounding  those  of  the 
empirical  spectra  are  used  to  compute  binomial-weight  MTI  improvement  factors, 
under  the  assumption  of  equal  -3  dB  widths  for  all  spectra.  It  is  shown  that 
great  differences  in  improvement  factor  can  be  expected  under  the  various 
shape  assumptions.  The  physical  reasonableness  of  the  illustrative  spectral 
functions  is  assessed. 
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A  major  conclusion  is  that  because  of  the  historical  bias  toward 
a  Gaussian  clutter  spectrum  shape,  the  interpretation  of  clutter  spectral 
width  constants  may  have  been  in  error  by  as  much  as  a  factor  of  2.5. 

1.3  Discussion 


The  existence  of  clutter  spectra  having  peaked  shapes  near  zero 
modulation  frequency,  rather  than  rounded,  brings  into  question  the  accuracy 
of  all  previous  estimates  of  clutter  spectrum  widths  which  were  made  under 
the  Gaussian  assumption  for  shape.  Direct  measurement  of  clutter  spectrum 
bandwidth  with  filters  must  be  corrected  differently  for  filter  bandwidth 
under  Gaussian  and  non-Gaussian  assumptions.  Indirect  measurement  of  bandwidth, 
via  a  decorrelation  time,  not  only  requires  appropriate  account  of  shape 
effects  but  also  of  the  effect  of  noise  near  zero  lag.  It  seems  probable 
that  some  pervious  measurements  reported  for  spectral  widths  may  be  in  error 
by  as  much  as  a  factor  of  2.5,  depending  on  how  the  estimate  of  width  was 
made.  This  error  is  somewhat  larger  than  one  would  expect  spectra  to  vary 
under  the  same  nominal  observing  conditions.  Without  a  thorough  review  of 
the  procedures  used  in  previously  reported  clutter  bandwidth  surveys  to 
transform  primary  measured  data  to  the  standardized  width  constants,  much 
more  about  the  bias  of  error  cannot  be  said  at  this  time. 

In  the  theoretical  modeling  of  clutter  spectra  presented  herein, 
some  variation  of  wind  shear  correlation  was  investigated  to  see  what  its 
effect  was  on  the  tails  of  modulation  spectra.  Over  the  range  of  correlation 
parameter  values  allowed  by  observed  wind  field  data,  the  predicted  spectrum 
tails  still  fall  off  faster  than  observed  clutter.  It  is  interpreted  that  a 
scatterer  lifetime  model  is  needed  in  addition  to  the  motion  model  to  add 
power  at  the  frequencies  to  the  theoretical  total  spectrum,  along  the  lines 
suggested  by  the  signal  analyses  in  Reference  2. 

The  existence  of  high-frequency  tails  in  the  spectra,  compared  to 
the  Gaussian  prototype,  has  implications  for  radar  Mil  design  and  performance 
prediction.  In  Appendix  B  are  shown  results  of  calculations  of  improvement 
factor  for  )-pulse  and  4-pulse  binomial  weighted  MTI  processors  under  the 
assumptions  of  clutter  spectral  shape  from  Gaussian  to  the  extreme  f**  (-2) 
and  the  same  -3  dB  bandwidth  for  each  analytic  function.  There  it  can  be 
seen  that  the  effect  of  the  tails  is  to  limit  the  performance  achievable 
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with  U  pulses  to  about  that  with  3  pulses,  and  to  temper  expectations 
about  the  improvement  achievable  to  lower  levels  than  the  ideal  Gaussian 
would  allow.  Given  such  hypothetical  non-Gaussian  spectra,  increased 
improvement  can  only  come  using  a  bank  of  bandpass  filters  appropriately 
weighted  to  the  spectrum  or  a  non-binomial  weighted  MTI. 

Appendix  C  examines  the  hypothetical  spectra,  each  of  which,  except 
the  Gaussian,  has  some  feature  relatable  to  the  empirical  data,  for  reason¬ 
ableness  in  regard  to  physical  considerations.  Finiteness  of  power  in  the 
process,  causality,  and  existence  of  derivatives  are  considered.  It  is  con¬ 
cluded  that,  of  the  set  examined,  only  a  spectrum  of  the  form  exp (-f **01)  , 
where  a<l,  passes  the  tests.  For  values  of  a  in  the  neighborhood  of  0.85 
that  function  also  approximately  represents  the  empirical  data. 


2.0 


EMPIRICAL  CLUTTER  SPECTRA 


Wide  dymanic-range  spectra  have  been  published  Cor  back-scattering 
from  foliage  and  from  sea  waves.  None  has  been  found  published  for  precipi¬ 
tation.  All  of  the  coherent  spectra  found  so  far  for  precipitation  have 
usable  dynamic  range  of  the  order  of  10  dB  and  uncertain  frequency  resolution. 
Here  are  presented  three  foliage  clutter  examples  and  one  sea  backscatter  and 
one  sea  f onward-scat  ter  example.  All  of  the  spectrum  plots  are  of  power  per 
unit  bandwidth  in  dB  against  modulation  frequency  on  a  linear  scale. 

2.1  Land  Clutter  Spectra 

Figure  1  shows  two  coherent  spectra  of  foliage  return  at  X  band 
[2],  The  curves  were  normalized  in  frequency  to  cross  at  their  -15  d3  points. 
The  publishing  authors  had  plotted  the  spectra  against  the  logarithim  of 
frequency,  thus  observing  the  simple  shape  and  especially  the  slope  near  zero 
frequency.  Wind  speeds  were  light,  5  kt  or  less  in  all  three  cases.  The 
cospectra  and  cross-spectra  were  computed  separately  and  showed  that  Doppler 
modulation  dominated  at  levels  above  the  -15  dB  point  but  amplitude  modulation 
contributed  below  that  level.  Amplitude  modulation  is  presumed  to  be  caused 
by  rotation  of  the  scatterers  in  the  wind. 


Another  example  of  radar  baokscatter  from  foliage  is  shown  in 
Figure  2  in  the  form  originally  published  [3] .  In  this  case,  the  frequency 
analysis  of  non-coherent  video  was  done  with  a  rectangular  time  window, 

which  produces  sidclobes  in  frequency  down  only  13  dB  peak  near  the  mainlobe 
and  of  the  order  of  -42  dB  at  f-40  Hz  and  -48  dB  at  80  Hz.  The  spectrum 
from  Figure  2,  which  is  plot  tea  against  fre«ucncy  on  a  logarithmic  scale, 
is  replotted  on  a  linear  frequency  scale,  in  Figure  3. 

A  second  example  of  noncoherent  spectrum  of  land  return  is  show 
in  Figure  4  and  5,  the  first  against  log  frequency  and  the  second  against 
linear  frequency  [4] ,  Again  there  is  a  tailing  departure  from  simple  expon¬ 
ential  but  at  a  much  lower  level  chan  predicted  in  Appendix  A,  and  the  slope 
near  zero  frequency  does  not  level  in  an  obvious  manner. 

2 . 2  Sea  Clutter  Spect  ra 

Figures  0  and  7  show  two  examples  of  non-coherent  spectra  of  radar- 
sea  return  (5).  Here,  over  the  somewhat  limited  dynamic  range  of  15  to  18 
dB,  the  data  can  bo  viewed  as  definitely  non-Gauss inn.  But  whether  the  tail 
of  the  spectrum  at  high  frequencies  falls  off  slower  as  suggested  by  the  other 
empirical  cases  cannot  be  known. 

The  remaining  empirical  spectrum  example  is  shown  in  Figures  8  and 
9.  There,  the  coherent  spectrum  of  L-band  waves  forward-scattered  from  a 
rough  water  surface  is  shown,  along  with  the  shape  of  the  analysis  filter  (6). 
Figure  8  details  the  spectrum  over  the  upper  30  dB,  and  Figure  9  shows  it  i'or 
a  wider  dynamic  range.  The  broadening  of  the  spectrum  below  -90  dll  is  believed 
caused  by  a  combination  of  interference  and  intermodulation  in  the  input  of  the 
spectrum  analyzer.  The  suggestion  is  strong  here  of  an  exponential  shape 
down  50  dB  or  so.  followed  by  a  departure  from  the  constant  slope  toward  higher 
power . 
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MODULATION  FREQUENCY  (CYCLES  PER  SECOND) 
Noncoherent  spectrum  of  radar  sea  backscatter  [5]. 


MODULATION  FREQUENCY  IN  HZ 


Coherent  spectrum  of  L-band  wave  reflected  from 
rough  water  surface  at  vertical  incidence  [6]  . 


MODULATION  FREQUENCY  IN  HZ 


igure  9,  Coherent  spectrum  of  L-band  wave  reflected  from  rough 


water  surface  at  vertical  incidence  [6] 
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A  MOTION  MODEL  FOR  CLUTTER  SPECTRA 
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Clutter  scatterer  velocities  seem  to  be  Gaussian  distributed. 

If  their  motion  were  independent,  one  would  expect  that  the  spectrum 
of  backscatter  signals  would  also  be  Gaussian.  This  logic  may  underlie 
the  historic  interpretation  of  observed  spectra  as  being  Gaussian  shaped. 

But,  there  is  strong  evidence  that  motion  independence  does  not  exist.  Where 
the  correlation  of  scatterer  motions  has  been  measured,  it  has  been  found 
that  the  scatterer  velocities  are  approximately  Wiener  distributed. 

Here  we  define  a  Wiener-distributed  velocity  process,  establish 
an  empirical  basis  for  its  reasonableness  in  one  case  (a  wind  field) ,  and 
compute  coherent  Doppler  specta. 

3.1  The  Wiener  Process 


A  random  process  x(t)  is  said  to  be  a  Wiener  process  if  it  satis¬ 
fies  these  conditions  [7] : 

(a)  x(t) ,  t>0,  has  stationary  independent  increments  for 
incremental  t; 

(b)  For  every  t>0,  x(t)  is  normally  distributed; 

(c)  For  all  t>0,  E(x(t))  =  0; 

(d)  x(0)  =  0. 

Var  [x(t)  -  x ( s ) ]  =  Oz  |t-s|,  (3) 

in  which  the  simple  parameter  a2,  plus  the  four  properties  above,  completely 
describe  the  process. 
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3.2 


An  Example  of  Scatterer  Motion 


Atmospheric  winds  containing  precipitation  seem  to  have  properties 
close  to  the  Wiener  process.  Consider  these  definitions. 


•  Air  moves  horizontally  with  an  average  vector  u  oriented  along 

the  x  axis.  Thus,  u  =  "iT  x 

x 


•  Turbulence  causes  the  local  wind  components  to  vary,  in 
particular  with  displacement  of  observing  position  along 
the  z(vertical)  axis.  Thus,  u  =  u^(z)x  +  u^(z)  y 


•  The  longitudinal  wind  speed,  u^  (z),  has  both  constant  and 
fluctuation  components.  Thus,  u^  (z)  =  u^_  (z)  +  v^(z) . 


In  the  transverse  direction,  u  (z) 

y 


Vy  (z) 


•  A  Structure  Function  is  the  mean  of  the  prod  ;et  of  spatial 
difference  of  two  wind  components.  Thus,  for  the  longitudinal 
co-structure  function  for  lag  is  the  vertical  direction, 


D 

xx 


OSl) 


Also,  D  (6z) 

yy 


^ux(z+6z)  -  ux(z) | 


£uz(z+  <5z)  -  Uy(z)  J2 


(4) 

(5) 


and  the  general  case  is 


V  )  =  h("+iV  "  ui(vl  '  uj  riv] 


•  Wind  shear  is  the  difference  in  a  component  of  the  wind  caused 
by  displacement  normal  to  the  component's  direction.  Thus, 
"vertical  wind  shear"  is  the  term  denoting  the  vector  cluinge 
in  horizontal  wind  for  displacement  along  the  horizontal 

— S  -N  -» 

axis:  Au(z)  =  Au  (z)x  +  Au  (z)y. 

_ x  y 

Also  Au^(z)  =  Aux(z)  +  Av^(z)  and  Auy(z)  =  Avv(z). 
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Observers  of  the  atmosphere  have  established  these  properties 
of  wind  shear.  Essenwanger  [8]  has  established  that 

(a)  the  fluctuation  components  of  the  wind  (turbulence) 

are  Gaussian,  with  the  standard  deviation  of  the  longi¬ 
tudinal  component  about  2.5  times  that  of  the  transverse 
component . 

(b)  The  fluctuation  components  of  the  vertical  shear  are 
Gaussian,  with  equal  standard  deviations. 

(c)  The  average  vertical  shear  is  about  one  third  the  standard 
deviation  of  the  longitudinal  component . 

(d)  The  standard  deviation  of  vertical  shear  increases  pro¬ 
portional  to  the  square  root  of  vertical  separation. 

One  of  Essenwanger ' s  data  sets  is  plotted  in  Figure  10.  There  the  logarithm 
of  shear  variance  is  plotted  against  the  logarithm  of  vertical  separation. 
The  slope  is  seen  to  be  about  1.1  for  separations  less  than  400m,  about 
0.7  above  400m,  with  an  average  of  1.0  over  the  maximum  1000m  span  plotted. 

Melnichuk  [9]  analyzes  shear  in  terms  of  the  longitudinal  and 

transverse  structure  functions.  He  finds  for  the  longitudinal  component 
near  the  ground  a  slope  of  1.2  below  a  break  at  400m  vertical  separation 

and  about  0.5  above,  up  to  his  maximum  analyzed  value  of  760m.  But  at  that 

low  altitude  (base  of  observations  was  130m  altitude) ,  the  transverse 
component's  slope  broke  at  150m  separation,  decreasing  from  1.0  below 
to  nearly  0  above. 
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Figure  10.  Variation  of  Vertical  Shear  Variance  with  Separation  [8] 


The  slope  referred  to  above  is  the  parameter  q  in  the  following 
empirical  laws: 

°vU)  =  A  Z<1 

(0) 

D (z)  «  C  2^ 

Both  authors  referred  to  above  heartily  subscribe  to  such  laws. 

3 . 3  Doppler  Spectrum  of  a  Wiener  Velocity  Process 

Let  a  radar  beam  illuminate  wind~driven  scatterers.  Choose 
the  center  of  the  beam  as  a  reference  point,  at  position  z-z^,  Assume 
that  the  wind  velocity  relative  to  the  average  wind  in  the  radial 
direction  is  representable  by  a  Wiener  process  in  the  variable  z.  Then 
v(z)  -  v(z  )  is  normal  with  a  2cro  mean  and  a  variance 


Var[v(z)  -  v (zo>  ]  =  a2  Iz-zJ.  (7) 

The  parameter  o  is  the  single  one  necessary  to  define  the  distribution 
of  this  shear  process. 

Assume  further  that  Che  average  reflectivity  of  scatterers  in 
a  layer  (a,  z+dz)  is  constant  with  z.  Then,  the  clutter  power  moving 
with  doppler  velocity  (v,  v  +  dv)  is 

S(v)  dv  »  dv  j  C(v,  z)y  G(z)  dz ,  (8) 


in  which y is  the  reflectivity  in  height  interval  dz  and  f(v,z)  is  the 
probability  density  function  of  v,  given  z.  G(z)  is  the  radar  beam 
shape  function. 
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The  pdl'  of  v  is 


» 


1 

f  (v)  «  -  exp 

V  /2?r  o 

v 


(9) 


in  which  is  the  variance  of  v,  here  simplified  to  denote  the  difference 

in  speeds,  v  *  v(z)  -  v(z  ).  the  variance  o  2=  a2 I z  -  z  I  is  as  in  Formula 

0  V  O 


Let  the  beam  shape  he  Gaussian, 


G(z)  **  exp 


so  that  its  3  dB  point  occurs  at 


ZQ .  5 


v'T  In  2  o 


1.1774  a_. 
is 


(10) 


(11) 


and  its  fuli  width  at  half  power  Is  2*/ln  2  «*  2,355o  .  The  beam’s 

one-dimensional  area,  or  effective  full  width,  is 

W  =  /HI  U2) 

and  the  total  scattered  power  is 

RCS  =  W  Y  =  yoB  (13) 
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The  function  S*  (v)  =  10  log^g(S(v))  is  plotted  in  Figure  11  For  three  values 
of  a,  o*  *  0.25,  1  and  4  mVsecVm.  Also  plotted  there  is  a  Gaussian  function 
having  small-lag  decorrelation  and  total  power  equal  to  t.hoso  values  for  the 
S(v)  (solid)  curve  for  which  o=l. 

Two  features  of  the  function  S(v)  are  noteworthy.  First,  the  slope 
near  v“0  is  not  zero,  in  marked  contrast  to  the  Gaussian  function.  Second, 
the  slope  is  nearly  constant,  varying  only  2:1  over  40  UN  change  in  the  spectral 
density . 


An  approximate  expression  for  S*(v)  is: 
S*(v)  - 


.  f~v  -  0,3  (20r<j^)"2/3  v-1 

[/ ?  r(1"*> '%  -J-] 


+  10  log 


(15) 


This  expression  matches  Formula  14  numerically  within  0.3  dB  ever  the  top 
70  dB  of  the  spectrum.  Although  not  visible  in  Figure  2  t he  slope  at  v^O 
is  zero,  and  the  radius  of  curvature  is  quite  small.  So  by  the  time  v  reaches 
a  few  percent  of  tV  spectrum  width  at  half  peak  power  density  the  slope  is 
significantly  non- zero.  However,  except  very  near  v=0,  the  slope  of  the  spectrum 
is  of  the  order  of  -5  dll  per  v  unit  for  c=l,  (see  Formula  15). 


Spectral 

Density 

(dB) 


Relative  Doppler  Frequency 

Figure  11.  Spectrum  of  Doppler  modulation  by  a  Wiener-distributed 
velocity  field. 
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3.4 


Other  Shear  Variance  Functions 


Define  the  variance  law  for  the  wind  shear  to  be 

Var[v(z)  -  v(zq))  =  a2  5zq  (16) 

where  the  exponent  q  is  now  arbitrary  but  of  the  order  of  one,  and  <5z  * 

| z— zq  j .  Using  this  law  instead  of  Formula  7  and  evaluating  Formula  8 
for  the  spectrum  gives  curves  plotted  in  Figure  12.  values  of  the  exponent 
less  than  1  are  more  curved,  tending  toward  the  Gaussian  shape.  Values 
greater  than  1  tend  toward  exponential,  and  a  value  of  1.6  yields  a  nearly 
straight  exponential  spectrum. 


12.  Spectrum  of  Doppler  modulation  by  scatterers  in  a 
modified  Vier.er-like  velocity  field. 
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APPENDIX  A 


NON-COHERENT  SPECTRA  OF  EXPONENTIAL-SPECTRUM  CLUTTER 


A. 1  Square-Law  Envelope  Detector 


Assume  a  clutter  signal  with  a  coherent  spectrum  of  the  form: 


S(f) 


=  A  exp  <  - 


(Al) 


After  the  signal  is  square-law  "envelope  detected'',  the  spectrum  of  the  non¬ 
coherent  video  signal,  neglecting  the  d.c.  component,  is  [a-i! 


Sn(f ) 


—  CO 


(A2) 


Using  Figure  A-l  as  a  guide,  this  integral  breaks  into  3  parts. 


ds 


ds 


ds 


(A3) 


A-l 


for  the  case  f>0. 


ds 


exp  t-  2s/f1>  ds 


s  (f)  =  -  A2  exp  (-f/f,)  /  exp  {+2s/f , }  ds 

n  7T  1  J  ■L 

— ca 

c 

+  i  A2  exp  (-f/f j 

QJ  e 

+  •”  A2  exp  (f/f^)  J 

£  ( 

S  (f)  =  -  A2  exp  (-f/f,)  (f i/2)  [exp  {2s/fJ 

n  it  i  x  J 

•00 

+  —  A2  exp  (-f  / f , )  £ 

rr  i 

c 

+  i  A2  exp  (f/fjK-fj/2)  [exp  {-23/fj}  j 


S  (f)  =  f,A2  exp  {-f/f,}  /2ir 

+  fA2exp  {-f/f^}  I' t 
+  f  1A2exp{-f /f exp{-2f/f1)  /2tt 


S  (f)  =  "  A?'exp{-f/f  }  (f+f  ) 
n  tt  II 

2 

sn(f)  =  (l+f/fj^)  exp  {-f/fj^} 


(A4) 


A- 2 


r 

i 

i 

E  Formula  A4  is  plotted  in  Figure  A-2  along  with  the  coherent  spectrum 

(A-l)  for  the  case  =  1,  but  the  densities  have  been  normalized  by 
displacing  the  non-coherent  curve  upward  by  lOlog^Q(ir). 


Figure  A-l 


Integration  Regions  for  Formula  A3. 


A-3 


Figure  A-2.  Noncoherent  Spectrum  (Square-law  detector) 
of  Exponential-Spectrum  Clutter. 
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A- 4 


Frequency 


LINEAR  ENVELOPE  DETECTOR 


For  che  linear  envelope  decector,  no  formula  for  Che 

video  spectrum  exists  directly  in  terms  of  che  coherent  spectrum.  Rather, 

the  auto  correlation  functions  muse  be  used.  Papoulis  derives  a  formula 

for  che  autocorrelation  of  the  linear-detected  video  signal  R^(t)  in  terms 

of  the  autocorrelation  function  of  the  coherent  signal  R  (t)  when  that 

c 

signal  is  a  Gaussian  distributed  stationary  process  [A-l], 


sin(a)  -  R  (t)/R  (0) 
c  c 


(A-5) 


R^Ct)  *  =  Rc(0)  fcos(a)  +  a  sin(a)) 


Note  that  because  the  coherent  signal  has  zero  mean,  as  T  varies  from  zero 
to  large  values,  a  varies  from  tt/2  to  approach  zero.  Thus,  one  sees  that 
R,  (T )  contains  a  dc  component  such  that  for  large  T,  R^(t)  approaches  2/tt 
Rc(0) .  3y  subtracting  this  term  from  the  autocorrelation,  the  auto  co- 
variance  of  che  video  spectrum  is  found. 


CL(x)  *  |  Rc(0)  (cos (a)  +  a  sin (a)  -  l) 


(A-7) 


Consider  a  coherent  signal  whose  spectrum  is  exponential.  Then  it  has  an 
autocorrelation  function  of  the  form: 


R  (T)  -  R  (0). 


1  +  02tV2 


(A-8) 


in  which  che  width  constant  C  has  been  so  defined  by  analogy  with  c'ne 
classical  Gaussian  parameter.  The  Fourier  transform  of  A-3,  which  represents 
che  coherenc  specerum,  is  [e.g.  A-3] 


(O  -  R  (0) 


7 rr 


exp  (2tt  y/T  f/*). 


A-5 


No  formula  exists  for  the  Fourier  transform  of  Formula  A-6  or  A- 7 
under  a  substituation  of  Formula  A-8  in  A-5.  Therefore,  the  spectrum  of  the 
non-coherent  video  signal  must  be  evaluated  numerically.  This  has  been  done 
for  both  the  coherent  and  non-coherent  aueocovariance  functions,  and  the  results 
are  plotted  in  Figure  A-3.  A  value  of  T^l  was  used  in  Formulas  A-8,  and 
the  frequency  scale  in  Figure  A-3  corresponds  to  this  choice. 

A-l  The  Fourier  Integral  and  Its  Applications,  A.  Papoulis, 

McGraw-Hill  Book  Company,  Inc.,  1962 

A- 2  Probability,  Random  Variables,  and  Stochastic  Processes , 

A.  Papoulis,  McGraw-Hill  Book  Company,  Inc.,  1965 
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APPENDIX  B 


NT I  IMPROVEMENT  FACTORS 

Ward  quotes  [  B-l]  the  IEEE  definition  of  MTI  Improvement  Factor 
(IF):  "The  signal- to-noise  ratio  at  the  output  of  the  system  divided  by 
the  signal-to-r.oise  ratio  at  the  input  of  the  system  averaged  uniformly  over 
all  target  radial  velocities  of  interst."  Here  binomial-weighted  MTI  cancellers 
are  assumed  to  constitute  the  "system".  The  IF  for  cancellers  operating  on 
3  and  A  pulses  is  given  by  the  following  expressions.  [B-2] 

3- Pulse  Canceller: 

IF  -  (6  -  8  p(T.)  +  2  ,'(2Tr))/6  (B-l) 

4- Pulse  Canceller: 


IF  *  (20  -  30  o(Tr)  +  12  p(2Tr) 


C  )  I -w 


(B-2) 


The  IF  for  these  two  filters  will  be  compared  when  fed  with  four- 
different  hypothetical  clutter  spectra.  The  typifying  spectra  chosen  are 
listed  in  Table  B-l.  They  have  all  been  defined  to  have  the  same  -3  dB  point. 
The  shape  of  the  exponential  spectrum  near  zero  modulation  frequency  is  typical 
of  empirical  spectra.  The  tails  of  empirical  spectra  at  higher  modulation 
frequencies  lie  between  the  f  ’  and  f  ''  examples.  The  Gaussian  spectrum  is 
included  because  of  its  historical  importance. 


The  results  of  computing  IF  using  Formulas  B-l  and  B-2  are  shown  in 
Figures  B-i  and  B-2,  respectively.  The  abscissa  variable  is  the  product  of 
interpulse  period,  Tf,  and  the  parameter  ,  the  -3  dB  point  of  the  spectrum. 
By  comparing  the  corresponding  curves  for  the  3-pulse  canceller  with  those 
for  the  4-pulse  case,  these  features  are  found.  For  the  Gaussian  spectrum, 
one  sees  the  distinct  difference  in  slope  of  IF  vs  T^  ^  beginning  at  about 
Tru.j  -  1.  For  this  spectrum,  large  difference  in  IF  would  be  seen  from 
cancellers  with  3  and  4  pulses,  when  is  small. 


B-l 


For  Che  exponential  spectrum,  there  is  also  a  change  in  slope, 
but  it  begins  only  below  about  T^io^  =0.2.  Because  of  the  distinctly 
different  shape  of  the  peaks  of  the  spectral  functions ,  the  IF  curve  for 
the  exponential  spectrum  is  shifted  toward  lower  values  of  by  a  factor 

of  2-2. 


Both  the  f  2  and  f  "  curve  show  no  change  in  slope  with  increase 
in  number  of  pulses  cancelled.  These  spectra  bound  the  expected  shape  of 
tails  of  empirical  clutter  spectva.  If  the  optimistic  f  tail  is  assumed, 
i.  binomial-weighted  4-pulse  MTI  canceller  increases  IF  by  about  3  dB  for  all 
below  about  0.1,  compared  to  a  3-puise  canceller. 

Two  conclusions  appear  warranted  from  this  information. 

•  Misinterpretation  of  the  shape  of  clutter  spectra  peaks  in 
the  past  may  have  caused  both  an  incorrect  estimate  of  the 
appropriate  spectrum  standard  deviation  constant  for  use 
Ln  MTI  calculations  and  an  optimistic  assessment  of  MTT  IF 
performance  that  would  be  aoheived  in  given  clutter  conditions. 

♦  Binomial  weights  for  MTI  cancellers  are  a  poor  choice  to  deal 
with  spectva  having  enhanced  tails  at  higher  modulation 
frequencies . 
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'ABLE  B-I  FOUR  EXEMPLARY  SPECTRUM  FUNCTIONS 


Improvement  Factor  for  a  3-pulse  Binomial-Weight  Canceller 
Spectrum  Types  are  Defined  in  Table  B-2 


APPENDIX  C 


PHYSICAL  ACCEPTABILITY  OF  HYPOTHETICAL  SPECTRUM  MODELS 

It  is  convenient  to  represent  empirical  spectra  in  terms  of  approxi¬ 
mating  functions  that  have  nice  mathematical  properties.  When  talking  about 
random  physical  processes  and  representing  their  fluctuation  behavior  in  terms 
of  these  selected  analytic  spectral  functions,  one  should  take  care  that  those 
functions  do  not  introduce  future  difficulties,  or  at  least  be  aware  of  the 
difficulties.  Examples  of  difficulties  include:  infinite  power  in  the  process; 
causal  contradiction;  and  nonexistence  of  derivatives  of  the  process.  Look 
at  each  of  these  in  turn. 

We  want  to  watch  out  for  implication  of  infinite  power  in  a  process. 
Empirical  spectra  of  certain  noise  processes  have  appeared  to  have  a  form 

S(f )  =  so/|f|.  (C-l) 

But,  this  spectrum  has  no  Fourier  transform  because  of  convergence  difficulty 
near  the  f=0  limit.  It  is  implied  that  the  power  from  a  source  with  such  a 
spectrum  is  infinite,  or  that  the  spectrum  fails  to  apply  at  some  small  value 
of  f.  When  attempts  to  find  the  cutoff  frequency  appropriate  for  some  processes 
have  been  made,  experimenters'  patience  was  exhausted  before  the  -  cutoff  was 
found.  [C-l]  Another  example  of  the  contradiction  about  power  is  the  white- 
noise  spectrum,  S (f )®coustant,  0<f<°c.  This  leads  to  the  classic  ultraviolet 
catastrophe,  [c-2]  Clearly  the  solution  to  both  of  these  dilemmas  is  to 
truncate  the  spectrum. 


The  notion  of  causality  implies  the  following.  If  the  input  of  a 
linear  system  has  been  zero  up  to  some  time  t,  then  its  output  has  also  been 
zero  for  all  time  up  to  t.  If  the  power  spectrum  of  the  output  of  the  system 
is  S ( f )  when  its  input  is  broadband  "oise ,  then  a  necessary  and  sufficient  condi¬ 
tion  that  the  system  is  causal  is 


C-l 


(C— 2) 


U&  df  <- 

1  +  f2 

(This  is  a  modification  of  the  Paley-Wiener  criterion,  applied  here  to  the 
power  spectrum  rather  chan  to  the  amplitude  spectrum.  Jc-3,-4,-5] ) .  There 
is  certainly  no  argument  about  the  need  for  causality  in  clutter  signals.  At 
some  prior  time  the  radar  was  not  transmitting.  (However,  the  clutter  as  a 
mechanical  process  may  have  been  going  for  a  longer  time.)  But  no  argument 
attempting  to  extend  time  works  when  we  recognize  that  clutter  is  not  stationary, 
and  its  statistics  behave  approximately  in  a  stationary  way  only  over  restricted 
sample  intervals.  One  might  well  argue  about  whether  an  empirical  random  process 
such  as  radar  clutter  need  be  representable  by  broadband  noise  passed  through  a 
linear  network.  Again  here  the  argument  depends  on  our  need  for  believing  in 
stationary  properties.  If  we  insist  on  using  the  notion  of  spectrum  to  represent 
an  expectation,  rather  than  a  deterministic  sample  of  frequency  content,  then 
we  must  accept  the  consequencies  of  the  condition  of  wide-sense  stationarity . 

One  of  these  is  that  a  network  exists  which,  when  fed  with  broadband  noise, 

has  as  its  output  a  proper  realization  of  the  process  whose  spectrum  is  the  one  in 

Question.  ^C— 3^J .  Thus,  the  condition  in  Formula  C-2  should  apply. 

Consider  the  implication  of  the  nth  derivative  of  a  clutter  signal 
being  zero.  Then  all  of  the  n-k  derivatives,  k=l..n,  are  deterministic,  not 
random.  Therefore,  for  the  clutter  signals  to  be  random  processes,  all  of  their 
derivatives  must  exist.  The  autocorrelation  function  of  the  derivative  of  a 
random  process  X  is 


and  the  power  of  the  derivative  is  Therefore,  we  require  that  the 

even-order  derivatives  of  the  autocorrelation  function  at  t=0  all  exist. 
Furthermore,  if  the  derivative  is  to  have  finite  power. 


f  f2n  S(f)df  <  °°  (C-4) 

-  CO 

must  apply  for  all  n.£c=3j 

Now,  consider  simple  spectral  types,  including  those  used  in  Appendix 
B  as  models  for  MTI  calculations,  and  apply  the  3  tests  above  to  them.  Table 
C— 1  summarizes  the  results. 

We  learn  t  i  tolerate  the  white-noise  model  because  our  attention  is 
always  centered  on  the  output  of  some  bandlimiting  system,  and  the  system  itself 
tells  us  where  and  how  to  truncate,  or  roll  off,  the  spectrum  and  thus  relieve 
the  contradiction  about  power.  Representation  of  empirical  spectra  with  a  modi¬ 
fied  white  model  always  includes  the  fact  of  a  rolloff,  if  not  the  correct 
shape  of  it.  The  1/f  spectrum  is  more  difficult  to  work  with  for  two  reasons. 
The  lower  limit  at  which  it  would  be  physically  reasonable  to  truncate  the 
spectrum  is  generally  obscure.  And,  the  variance  of  the  process,  and  hence 
the  power  ascribed  to  it,  is  proportional  to  the  measurement  time.  Thus,  such 
processes  are  not  stationary  and  the  notion  of  spectrum  as  expectation  is  hard 
to  live  with. 

The  Gaussian  spectrum  is  far  from  satisfying  the  causal  test,  but 
the  Exponential  is  close.  That  is,  spectra  of  the  form 

S(f)  =  Sq  exp(-fa)  (C-5) 

are  acceptable  so  long  as  a<l.  Note  that  when  a<l,  sucli  spectra  have 
high-frequency  tails  qualitatively  in  agreement  with  empirical  clutter 
spectra,  in  that  both  deviate  from  exponential  (ct=l)  in  the  same  direction. 


\ 


C-3 


The  last  two  spectral  types  in  Table  C-l  have  difficulty  in  meeting 
the  derivative  test.  The  derivatives  have  infinite  power  above  n=l  and  2, 
respectively. 

Thus,  only  one  type  of  spectrum  discussed  so  far  satisfies  all 
three  tests,  the  modified  exponential  defined  in  Formula  C-5. 
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APPENDIX  D 


SELECTION  OF  SPECTRUM  WIDTH  PARAMETER 


For  the  case  of  precipitation  clutter,  the  vertical  shear 
properties  are  well  documented.  But  radar  designers  need  clutter  spectral 
properties  to  guide  selection  of  MTI  and  waveform  designs.  The  physical 
model  investigated  in  Section  3.3  provides  the  means  for  relating  the 
wind  shear  to  the  clutter  Doppler  spectrum  width.  But  the  path  needs 
the  clarification  that  follows. 

We  accept  Chat  for  some  stated  wind  conditions,  a  value  (or 
range  of  values)  for  the  vertical  shear  variance  can  be  agreed  upon. 

Then  by  Formula  7 , 

a v~(5z)  »o'l<5zl  (D-i) 

is  the  variance  of  the  velocity  difference,  v,  between  scatterers  separated 
by  distance  «5z,  in  terms  of  the  Wiener  parameter,  a.  Therefore,  it 
follows  that 

a2  =  av2  (5  z)  / 1  <5z  I  ••  (D-2) 

can  be  used  to  estimate  appropriate  values  of  the  parameter  a.  Alternatively, 
one  might  want  to  use  a  variation  based  also  on  Formula  6.  The  spectrum 
computation  in  Section  3.3  which  guides  the  next  steps  would  apply  correctly 
to  both  cases  provided  the  spectrum  slope  parameter  as  defined  by  Formula 
15,  is  interpreted  in  light  of  the  variations  seen  in  Figure  12. 


Suppose  that  an  exponential.- Like  spectrum,  such  as  the  ones 
plotted  in  Figures  11  and  12  and  which  are  associated  with  Wiener-like 
velocity  fields,  suit  the  radar  analyst's  needs.  By  Formula  15,  the 
spectrum  in  the  vicinity  of  v  =  0  is  for  =  1 


S*(v) 


-  —  v  +  constant . 

a 


(D-3) 


\  j  1 

If  v  has  units  of  meters/sec  and  a  units  of  m  /sec,  then  the  slope 

3/2 

constant  must  have  units  of  dB/(m  /sec)  in  order  that  S*(v)  nave  uni 
of  dB/(m/sec).  But  v  requires  calibration,  because  of  the  restrictive 
assumption  that  Che  beam  height  parameter  have  value  1.  Inspection 
of  Formula  14  shows  that  the  proper  calibration  is  made  by  replacing  v 
by  v/oQ.  This  may  be  done  for  Che  abscissas  of  Figures  11  and  12  and 
for  Formulas  15  and  D-3  above.  Thus, 


ts 


S*(v) 


- —  +  c(o,oB) 

OVOg 


(D-4 ) 


has  the  needed  scale  factor.  The  "hidden"  units  now  associated  with  the 

-  ,  1/2 
factor  5  are  dB/m 


at 


From  D-4  it  follows  that  the  3  dB  point  on  the  spectrum  occurs 


Thus,  by  Formula  D-2  it  follows  that 

v3dB  , 

and  in  Doppler  angular  frequency  space,  where 

U)  =  4tt  v/\  , 

&).,  =  2.4  ::/aJ  0  (5z)/ GW 6?) . 

3  B  v 


(D-5 ) 


(D-6) 


(D-7) 

(D-8) 
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Suppose  that  a  radar  has  a  3  dB  elevation  beamwidth  (for  2-way 
propagation)  of  t>3  .  Then  by  l-’ormula  li  (and  the  comment  following  it), 

a0  =  03/(21n(2))  =  0^/2.355  .  (D-9) 

At  a  radar  range  of  R,  the  height  of  the  beam  at  its  -3  dB  (2-way) 
points  is 


h  ^  =  R03  , 


from  which  it  follows  that 


=  ah  =  R°0  =  RV2’335  • 


By  combining  D-7  and  D-10,  the  formula 


to. 


2 . 4tt av(6z) 


/ 2 . 355’  X 


<03  = 


4 . 9/R0?  a  (Sz) 
_ 3  v 

X  /&' 


(D-10) 


(D-ll) 


(D-12) 


CD-13) 


Consider  the  following  numerical  example.  A  typical  value 
for  atmospheric  vertical  shear  standard  deviation  is  quoted  as 
5  m/sec/km^2  { D—  1  ] .  Thus,  0  ^  ( <5  z )  f/  tS  z1  =  5//10001  ~  0.158  m^2/sec.  For 
a  3.2-cm  radar  having  2-way  elevation-plane  beamwidth  of  2  degrees  viewing 
a  target  cell  at  50  Itm  range,  the  Doppler  bandwidth  for  widespread 
precipitation  is 


-  9.9  / 50 , 000  , 0349 1  ,158 

3  "  ,032 


=  1010  radians/sec. 


Rc  ference 

[D-l  1  S.  I,.  Valley,  Handbook  of  Geophysics  and  Space  Environments, 
Air  Force  Cambridge  Research  Laboratories,  1965. 
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